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SIGNIFICANCE  AND  EXPLANATION 


He  present  an  analysis  of  a  one-dimensional  model  of  single- junction 
semiconductor  devices  (pn- junctions  and  certain  resistors)  when  an  external 
voltage  is  applied  to  the  contacts.  The  model  has  the  font  of  a  system  of  six 
highly  nonlinear  first  order  ordinary  differential  equations  subject  to 
boundary  conditions  at  the  contacts  of  the  device.  The  system  is  singularly 
perturbed  (the  derivatives  of  some  components  are  multiplied  by  a  small 
constant/  the  so  called  singular  perturbation  parameter).  The  dependent 
variables  are  the  electrostatic  potential/  the  hole  and  electron  densities  and 
the  hole  and  electron  current  densities.  A  region  of  fast  variation  in  the 
electrostatic  potential  and  in  the  carrier  distributions  occurs  due  to  the 
singular  perturbation  character  of  the  problem.  This  region  is  in  the 
interior  of  the  device  (internal  layer)  and  represents  the  junction  between 
differently  doped  areas.  We  derive  formal  asymptotic  expansions  of  solutions 
as  the  singular  perturbation  parameter  tends  to  zero  and  we  prove  that  such  an 
expansion  'represents'  a  solution.  We  also  investigate  the  dependence  of  the 
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ABSTRACT 

In  this  paper  we  present  an  analysis  of  the  fundamental  one-dimensional 
semiconductor  equations  describing  potential,  carrier,  and  current  density 
distributions  in  single- junction  semiconductor  devices  when  an  external 
voltage  is  applied  to  the  contacts.  We  reformulate  the  model  equations  by 
appropriate  scaling  as  a  singularly  perturbed  two  point  boundary  value  problem 
for  a  system  of  nonlinear  ordinary  differential  equations.  The  right-hand 
side  of  the  system  has  a  jump  discontinuity  with  respect  to  the  independent 
variable  (space-coordinate)  representing  the  junction  between  differently 
doped  sides  of  the  device.  The  solution  components  are  assumed  to  be 
continuous  across  this  junction.  - 

/  v 

We  give  an  existence  proof  for  the  reduced  problem  (the  singular 
perturbation  parameter  is  set  to  zero).  The  discontinuity  of  the  right-hand 
side  of  the  system  produces  a  discontinuity  in  the  reduced  potential  and 
reduced  carrier  distributions.  This  creates  an  internal  layer  in  the 
corresponding  solution  components  of  the  singularly  perturbed  problem.  The 
current  distributions  have  no  internal  layer.  We  also  derive  the  (internal) 
layer  equations  and  give  an  existence  proof.  No  boundary  layers  occur. 

We  show  that  formal  expansions  actually  represent  (asymptotically) 
solutions  of  the  singularly  perturbed  problem  if  the  applied  voltage  is 
sufficiently  small,  and  we  investigate  the  dependence  of  the  total  current  on 
the  applied  voltage.  Numerical  computations  are  reported. 
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AN  ASYMPTOTIC  ANALYSIS  OF  SINGLE- JUNCTION 
SEMICONDUCTOR  DEVICES 


*  ** 

Peter  A.  Markowich  ,  C.  A.  Ringhofer  , 

***  *** 

E.  Langer  and  S.  Selberherr 

1.  Introduction 

In  this  paper  we  present  an  analysis  of  a  class  of  systems  of  ordinary  differential 
equations,  subject  to  boundary  conditions,  modelling  pn-junction  devices.  The  physical 
situation  is  as  follows.  A  semiconductor  (for  example  Silicon)  is  doped  with 
acceptor  atoms  (negative  ions)  in  the  left  side,  with  donor  atoms  (positive  ions)  in  the 
right  hand  side  and  a  bias  U  =  Uft  -  Uc  is  applied  to  the  contacts: 


cathode  contact 
Applied 

Cathode  potential 
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The  device  is  assumed  to  have  characteristic  length  2l(*  5  x  10  cm)  and  the  junction  is 
at  z  =  Z  e  (-Z,l)  (the  term  junction  refers  to  the  boundary  of  the  n  and  p  regions  as 
well  as  to  the  whole  device).  The  device  is  forward  biased  for  U  >  0  and  reverse  biased 
for  U  <  0.  The  physics  of  pn-junction  is  explained  in  Sze  (1969),  Ashcroft  and 
(1976)  and  R.  A.  Smith  (1978). 
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The  equations  describing  the  electrostatic  potential,  the  carrier  densities  and  the 


current  densities  within  the  device  (in  the  static,  one-dinensional  case)  are: 


(a)  ^  (n-p-INpts)  -  N^tz)))  Poisson's  equation 


(b)  n  “  —  ni  ♦  -  J  electron  current  relation 

D  qD  n  . 

n  n 


(1.1) 


(c)  p'  “  p<i'  ”  —g—  Jp  hole  current  relation 
P  q  P 


(d)  J'  -  qP(n,p,J  ,J  )  continuity  equation  for  electrons 
n  n  p 


(e)  J  m  qP(n,p,J  ,J  )  continuity  equation  for  holes 
P  n  P 


for  z  e  (-i, t)  ("'"  denotes  differentiation  with  respect  to  z)  subject  to  the  boundary 
conditions 

"i 

(a)  ♦(-<)»  UTtn  ■  (anode) 


(1.2) 


(b)  ♦(»)  -  U  in  ♦  U  (cathode) 

t  c 

(c)  n(±i)p(£t)  «  n| 


(d)  n(±t)  -  p(+i)  -  N*(it)  -  »T(+t)  • 

(See  Van  Roosbroeck  (1950).) 

The  dependent  variables  (with  units)  in  (1.1),  (1.2)  are 

♦  :  electrostatic  potential  (V) 

:  electrostatic  field  (Vcm  ') 
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n  :  electron  density  (cm  ) 
p  :  hole  density  (cm 

2 

J  :  electron  current  density  (A/cm  ) 
n 

2 

J  :  hole  current  density  (A/cm  )  . 

P 

All  parameters  in  (1.1),  (1.2)  (except  N  (z),  N*(z))  and  the  temperature  T  are  assumed 


to  be  constant.  Table  1  gives  the  physical  meaning  and  approximate  numerical  values  of 
these  parameters  at  T  *  300X  (room  temperature)  for  silicon. 


Table  ft  Parameters  for  Silicon  at  T  *  300k 


Parameter 

Physical  Meaning 

Numerical  Value 

'1 

elementary  charge 

10  ’’fts 

e 

permittivity  constant 

10-12AS/Vcm 

“n 

electron  mobility 

103cm2/Vs 

hole  mobility 

103cm2/Vs 

Dn 

electron  diffusion  constant 

25  cm2/s 

Dr 

hole  diffusion  constant 

25  cm2/S 

"i 

D  D 

intrinsic  number 

ID10™"3 

o  »  ~  — E 

n  p 

thermal  voltage 

0.025V 

is  the  density  of  electrically  active  acceptor  atoms  and  N*  Is  the  density  of 
electrically  active  donor  atoms  and 

(1.3)  C(z)  =  N*(z)  ”  ( z)  (cm  3) 

is  called  doping  (or  impurity)  profile.  For  the  pn-junction  C(z)  is  negative  for 
z  e  [-t,Z)  (p-side)  and  positive  for  z  S  (Z,tl  (n-side)  and  is  assumed  to  have  a  jump- 
discontinuity  at  z  ■  Z  (abrupt  junction).  He  also  investigate  the  less  important  and 
much  simpler  case  C(z)  >0  in  (-t.tl  (but  still  with  a  jump-discontinuity  at  z  •  Z). 
These  devices  are  called  n+n-  or  nn+- junctions  (depending  on  whether 
C(z, )  >  C(z2)  or  C(z,)  <  C(z2)  for  all  z^  «  [—  t,X) ,  z2  «■  (X,D). 

The  analysis  of  p+p  and  r>p+  junctions  (Ciz)  <0  on  [-1,  U  )  is  analogous  to  the 
analysis  of  n+n  and  nn+  junctions.  Only  n  and  p,  Jn  and  have  to  be 

interchanged  and  f  has  to  be  substi tilted  by  -fi. 

2  2 

The  scalar  function  R  e  C(  (0,«*)  x  r  )  in  (1.1)(e)  is  called  recombination  term,  it 
is  the  rate  at  which  electron-hole  carrier  pairs  are  generated  (R  <  0)  or  recombine 


! 

I 

I 

i 


I 


l 


i 


2 

(vanish)  (R  >  0),  R(n,p,0,0)  -  0  for  n,  p  such  that  p  •  holds  (equilibrium 
condition) • 

The  Shock ley-Read-Ha 11  ( SRH )  recombination  term 

2 


(1.4) 


_ nP  ~  _ 

RSRH("'P>  "  T  (n+n  )  r  (p+n.  ) 
p  i  n  i 


/  -3  -1, 

(cm  a  ) 
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describing  thermal  recombination,  where  t  ,t  { •>  10  s)  are  the  electron  and  hole 

n  p 

lifetimes,  is  widely  used.  Different  ways  of  modelling  R  (which  are  necessary  for  very 
large  |U|  )  are  given  in  Langer,  Selberherr  and  Mader  (1981)  and  SchSts,  Selberherr  and 
PBtzl  (1982). 


The  boundary  conditions  (1.2)  (c)  express  that  the  contacts  a  “  ±1  are  in  thermal 
equilibrium  and  (1.2)(d)  represents  vanishing  space  charge  at  the  contacts. 

We  only  admit  solutions  of  (1.1),  (1.2)  which  fulfill 

( 1 .5)  (a)  <i,n,p,J  ,J  6  c’d-Ml) 

n  p 

( 1.5)  (b)  e  c\ l-t,z)),  ♦’  e  c1  ( (z,tl ) 

(1.6)  n  >  0,  p  >  0  on  . 

(1.5)  comes  from  the  jump-discontinuity  of  C(z)  (♦H  cannot  be  continuous  if  n,  p  are 
continuous),  the  equation  (1.11(a)  has  to  be  fulfilled  for  the  right  hand  limit  and  for 
the  left  hand  limit  of  4i”  at  *  m  Z.  (1.6)  has  to  hold  on  physical  grounds  since  n,  p 
are  densitites . 

In  this  paper  we  scale  (1.1),  (1.2)  such  that  we  obtain  a  singular  perturbation 

X 

problem.  The  perturbation  parameter  (called  X  in  the  sequel)  is  equal  to  - —  where 
euT  1/2  1 
XD  “  (q  nuixIcYsYj)  ia  the  length. 

we  present  an  asymptotic  analysis  of  (1.1),  (1.2)  (for  X  small,  which  corresponds  to 
large  dopino  |C|).  The  discontinuity  of  C(z)  at  z  «  Z  produces  an  (internal)  layer  in 
the  fast  components  ^t^'intp.  Jn,  Jp  are  the  slow  components  (uniformly  c1  as 
X  ♦  0+). 
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We  derive  the  reduced  problem  (X  »  0),  the  layer  equations  and  give  existence 
theorems  for  both.  The  reduced  problem  has  the  form  of  a  two-point  boundary  value  problem 
with  interface  conditions  at  the  discontinuity.  Using  these  results  we  prove  an  existence 
result  for  the  full  problem  (1.1),  (1.2)  (for  sufficiently  large  doping  |C|,  that  means 
X  small)  assuming  that  the  recombination  rate  R  =  0  (corresponding  to  infinite  electron 
and  hole  lifetime)  and  that  |U|  is  small.  We  show  that  (for  X  sufficiently  snail) 
there  is  a  solution  of  (1.1),  (1.2)  whose  fast  components  are  close  to  the  sum  of  the 
(corresponding)  'reduced'  solution  components  and  the  layer  terms  and  whose  slow  components 
are  close  to  the  corresponding  'reduced'  solution  components.  Ho  layers  at  the  contacts 
occur  since  the  'reduced'  solution  fulfills  the  'reduced'  boundary  conditions. 

we  also  investigate  the  dependence  of  the  total  current  J  “  Jn  ♦  Jp  on  the  applied 
voltage  U  (J  is  a  constant  because  of  ( 1 . 1 ) (d) , (e) > ) .  It  turns  out  that  J  is 
asymptotically  (as  X  ♦  0+)  a  linear  function  of  U  if  C(z)  >0  on  (n+n  and 

nn+  junctions  are  resistors)  and  J  is  asymptotically  an  exponential  function  of  U  if 
C  changes  sign  at  Z. 

The  singular  perturbation  approach  to  pn-junction  modelling  was  suggested  by  many 
authors.  Vasilev'a  and  Butuzow  (1978),  Vasilev's  and  Stelmakh  (1977)  and  D.  Smith  (1980) 
investigated  a  much  simplified  model  (they  assume  that  the  current  densities  are  known 
instead  of  the  applied  voltage,  that  Z  ”  0  and  that  C(z)  is  odd)  and  prove  an  existence 
theorem  using  the  asymptotic  expansions.  The  authors  of  this  paper  analyzed  (1.1),  (1.2) 
(1982)  under  the  (pretty  unrealistic)  assumption  that  the  junction  Z  is  in  the  middle  of 
the  device  and  that  the  doping  profile  C  is  odd.  This  allows  to  reduce  the  internal 
layer  problem  to  a  boundary  layer  problem.  The  advantages  of  the  singular  perturbation 
approach  for  the  numerical  solution  of  (1.1),  (1.2)  is  also  explained  in  the  latter  paper. 

The  generalization  of  the  presented  theory  to  multilayer  structures  like  bulk-barrier 
diodes  (see  Langer,  Selberherr  and  Mader  (1981))  or  thyristors  (see  Sze  (1969))  is 
straightforward. 
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Thi«  paper  la  organized  as  follows.  In  Section  2  we  perform  the  scaling  and 
reformulation  of  <1.1),  (1.2)  as  a  singular  perturbation  problem,  in  Section  3  we  derive 
the  expansions,  prove  existence  theorems  for  the  reduced  problem  and  the  internal  layer 
problem  and  in  Section  4  we  give  the  existence  proof  for  F  =  0  the  full  singularly 
perturbed  problem  in  the  case  \  and  |0|  small.  Numerical  results  for  large  0  are 
demonstrated  in  Section  5. 


2.  Scalin' 


we  scale  the  dependent  variables  as  follows 


(2.1) 


*■ '  V  ‘ 


tJ 


tJ 


(2.2) 


S  0  qC  B  D  qC 
n 


where  C  *  max  |C(z)|  and  the  Independent  variable 
ze [-1,1] 


(2.3) 

Then  (1.1)  reads 


Z 

x  ■*  — 
1 


(a) 

>V 

•  n  -  p 

8 

a  8 

(b) 

n'  = 

nf  +  J 

8 

s  s  n 

(c) 

p:  - 

-p  if  -  J 

s  o'  m  « 

(2.4) 


,2  D  qC  0  oC 

(d)  j;  -  —  R(cna.  Cp  #  -V  *n  .  -V  ) 

a  D  c  s  s 

n 


,2  __  DgC  D  gC 

(e)  J*  «  -  -  r( Cn  , Cp  ,  J  •  j  ) 

P8  dc8"4  "a  1  ps 
P 


D(x)  =  (  |  D  ( x  )  |  <  1) 


for  -1  <  x  <  1  ("'"  denotes  now  differentiation  with  respect  to  x).  We  have  set 

(2.5) 
and 

(2.6) 


,2  XD  eUT 

*  7m7j  • 


The  boundary  conditions  (2.2)  are 

{a)  v-„  -  ♦  r 


2, 2  U 


s  T 

n  ( 1 )  u 


(2.f ) 


(b)  *,(1)  ■ 

y  \  T 


(c)  n  (+1 )p  (£1)  »  y4A4 
s  s 


(d)  n  (  + 1 )  -  p  (±D  *  D(  +  1) 


where 


2 

(2.9)  Y  = 

T 

holds. 

We  now  assume  that  the  recombination  term  R  is  such  that 

D  qC  D  qC 


(2.10) 


(cn  ,Cp  ,  -7 —  J  ,  ?  J  )  “  CS(n  ,p  ,J  ,J  »Y*> 

'  s  s  t  n  Jt  p  '  ssnp 


2  2  — 
holds,  where  S  e  C( (0,«)  x  R  x  (0,-))  is  independent  of  C. 

Dropping  the  index  s  the  problem  now  reads 


(2.11) 


(a) 

X2*" 

=  n-p-D(x 

(b) 

n'  = 

"**  +  Jn 

(c) 

P*  * 

-P*'  -  Jp 

(d) 

II 

-  c 

*-3 

S  (n,p,J 
n  n 

(e) 

r  •  = 

P 

-S  ( n , p, J 
P  ] 

-1  <  x  <  1 


with  S  =  —  S,  S  =  —  S  subject  to  the  boundary  conditions 
n  D  p  D  ,  ,  „ 

n  P  2.2  U 

(a)  *(-D  =  to(^Tn)  +  b7 


(2.12) 


(b)  +( 1 )  =  £n(Iy^)  +  ^ 

y  K  T 


(c)  n(±1  )p(±1  i  =  y4  X4 


(d)  n( ± 1 )  -  p(±1 )  =  D(±1 ! 


If  D  =  D  =  D  holds  we  have 
n  p 

(2.13) 


S  =  S 
n  P 


•  i.  a  u  * 

Under  this  assumption  and  x  =  x  =  T  we  get  for  the  SRH  term  with  P  =  “r 

n  p  .2 


,  4*4 
S  -  e  -  1  rtp-T  X 

n  p  0  ,  2,2  ' 

n+p+2y  X 


Generally,  the  equilibrium  condition  implies  that  { n,p, 0 , 0 , yX)  -  S^(n,p, 0, 0 , yX) 

4  4 

holds  for  n,  p  such  that  np  -  y  X  . 
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The  discontinuity  of  0  occurs  at 

(2.15)  X  -  j 

and  th«  conditions  on  the  solution  of  (2.11).  (2.12)  are 

(2.16)  (a)  *,n,p,J  ,J  e  c’(!-1,1]) 

n  p 

(2.16)  (b)  e  c’(  (-1.X) ),  e  c\  (x,u ) 

(2.17)  n  >  0,  p  >  0  on  (-1,1) 

(see  (1.5),  (1.6)). 

The  boundary  values  for  n,  p  and  41  can  be  computed  from  (2. 12)  (c) ,  (d) : 

(2.18)  n ( 1 )  -  ~  (nd)  +  /d(1)2  +  4-fV,  p(  1 )  -  ~  (-D(I)  ♦  /d(1)2  ♦  4y4A4) 

(2.19)  n(-1)  -  J  (0(-1)  +  /d(-1)2  +  4y4A4),  p(-1)  -  ^  (-D(-I)  ♦  /d(-1)2  ♦  4y4A4) 


(2.20) 


♦(-1 )  -  in 


2.2 

2r  A 


-D(-1) 


+  /d"('-1)2  +  4 


u 

♦  -A 


(2.21) 


♦  (1)  -  in 


-D(1)  *  /d(1)2  ♦  4y4A4 

L 


♦  ^  . 

°T 


For  X  small  the  problem  (2.11),  (2.18),  (2.19),  (2.20),  (2.21)  constitutes  a  singularly 

perturbed  two-point  boundary  value  problem. 

A  small  means  that  C  is  large  (assuming  that  i  is  constant).  In  practical  cases 
“17  2  —6 

C  >  10  such  that  X  <  0.4  *  10  holds.  For  the  following  analysis  we  assuste  for  the 
sake  of  simplicity  that  D(x)  is  Independent  of  X  (it  would  suffice  to  assume  that  D 
is  analytic  in  X),  that  means  that  the  doping  |C(z)|'  Increases  'as  a  whole'  as 
X  ♦  0  +  .  Actually,  the  asymptotic  analysis  presented  in  the  next  sections  requires  that 


(2.22) 


2.2 

and  since  y  X 


(2.23) 


min  )D(z)| 

ze [-i.il 

max  1 D( z ) | 
z«[-t,t] 


>>  A 


max  | C ( z ) | 

ze(-i,i) 


that 


n  <<  max  |C(z) I 

1  ze [-l.il 
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holds 


The  two  casss  ws  deal  with  now  ara 


D(x)  <  0  for  x  e  [-1,X];  D(x)  >  0  for  x  e  <X,1]  and 
I D( x) |  >  Da  for  x  e  [-1,1] 


which  corresponds  to  the  pn-junction  and 

(B)  D(x)  >  DB  >  0  for  x  «  1—1,11 

corresponding  to  an  nn+  or  n+n  junction-  In  both  cases  (A)  and  (B)  we  assume  that 

D(X+ ) ,  D(X-)  (we  use  the  notation  of  f(X±  “  11m  f(x)  in  the  sequel)  exist  and 

x*X± 

D(X+)  1  D(X-)  and  that  D  is  sufficiently  smooth  everywhere  else. 

The  analysis  of  the  scaled  problem  is  complicated  by  the  loqarithmic  blow-up  of  the 
boundary  data  of  il”  as  given  by  (2.20 ),  (2.21)  in  the  case  (A). 

The  potential  difference  of  the  contacts  is  given  by 


(2.24) 


<M-1>  -  *(1) 


o  ,  VX) 

u- 

T  T 


where  the  build-in-voltage  Ubi  (X )  (i.e.  the  voltage  due  to  doping)  is  given  by 


U.  (X)  t"  / - 5 - 

(2.25)  A-.tn  D(-11 

T  D  ( 1 )  +  /d(1)2+«Y4A4 


tn^0(iHP~i)T^  +  0,yV)'  d,-1)  4  0  (A) 


tnf~-p )  +  0(y4X4),  D(-1)  >  0  (B)  . 


U,  (X)  is  bounded  as  X  ♦  0+  in  the  case  (B).  Since  (2.11)  depends  only  on  ♦"  (and 

bi 

not  on  <|i)  we  can  therefore  remove  the  singularity  in  the  boundary  conditions  by 

substituting  ♦  by  -  tn  -  ■ — The  equations  (2.11)  remain  unchanged  and  the  new 

Y  a 

boundary  conditions  for  the  case  (B)  arei 


(2.26) 


*<-1)  -  »n[j  (D(-1)  *  /d(-1)2+4y4X4)]  +  ~ 

T 


(2.27) 


*<1>  -  tn(j  (D(  1 )  +  /d(1)2+4y4X4)]  +  — ■  . 

T 
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m.'  J 


3.  Expansions 

In  this  Section  we  apply  the  standard  approach  for  singularly  perturbed  boundary  value 
problems  to  the  semiconductor  problem.  We  assume  that  the  solution  has  a  formal 
asymptotic  expansion  in  X,  each  term  in  the  series  being  the  sum  of  a  uniformly  smooth 
function  and  layer  terms. 

A  problem  that  occurs  is  the  blow  up  of  the  boundary  values  of  $  in  the  case  (A), 
which  implies  that  'reduced'  boundary  conditions  (XpO)  cannot  be  defined  formally  for 
V-  For  the  derivation  of  the  expansions  we  set  in  the  case  (A) 

(3.1)  (a)  (a)  t(-1)  -  (b)  *<1)  - 

and  assume  that  if  ,  <>+  are  independent  of  X  (this  will  be  justified  later).  In  the 

case  (0)  'reduced’  boundary  conditions  for  (r  can  be  obtained  from  (2.26),  (2.27)  and  we 

sets 

UA  Uc 

(3.1)  (c)  *_  =*  In  D(-1)  ♦  jp,  (d)  =  In  D(1)  ^  ^  . 


We  sake  the  following  'ansatz'i 

(3.2)  (a)  t(x,X)  -  tlx)  +  *(0)  ♦  tt<*>  ♦  *r(4>  ♦  . 

(3.21(b)  n(x,X)  ~  n(x)  t  n(0)  ♦  n^tt)  ♦  nf(^)  ♦  ••••• 

(3.21(c)  p(x,X)  -  p(x)  ♦  p(g)  +  Pj(t)  ♦  pr  ( ♦)  +  ••••• 

(3.21(d)  J  (x,X)  ~  J  (x)  ♦  J  («)  ♦  J  <T)  ♦  J  (♦)  . . 

n  n  n  n .  n 

t  r 

( 3.2 )  (e)  J  (x.X)  **J  (x)  ♦  J  (e)  ♦  J  ID  *  3  (♦>  . 

P  P  P  Pt  Pr 

where  the  dots  stand  for  a  power  aeries  in  X  (starting  with  the  0(  X)-term)  whose 
coefficients  are  of  the  same  form  as  the  given  0(1)  terms.  The  fast  variables  are 


(3.3) 

(3.3) 

(3.3) 


(a )  o 


( b)  T  = 


x-X 

X 

x+1 


\  x  e  (-i,i) 


<c> 


♦  = 


-J 


I 


.  : 

tel 


u 
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The  functions  narked  with  1 - '  denote  the  reduced  solution  (of  order  zero),  ‘ " '  denotes 
the  internal  layer  terns  (at  x  ■  X)  (of  order  zero),  '«■'  denotes  the  layer  terms  (of 
order  zero)  decaying  from  the  left  boundary  x  «  -1  (with  index  t)  and  the  layer  terms 
(of  order  zero)  decaying  from  the  right  boundary  x  ”  +1  (with  index  r)  resp.  The 
boundary  condition 


(3.4) 

(a) 

♦(+“)  “  n<±») 

«  p(±*">  “  J  (+•)  “  J  (i») 
n  p 

=  0 

(3.4) 

(b) 

(-)  -  nt(») 

•  p, (»)  ”  J  (•)  “  J  (“)  = 
n  p 

0 

(3.4) 

(c) 

♦  (-“)•  n  { 

-  •)  -  p  (“  “)  “  J  (-«•)“ 

J  (-  -)  -  0 

r  r 

r  n 

p 

hold,  since  the  internal  layer  terms  are  regarded  as  functions  on  R,  the  left  layer  terms 
are  regarded  as  functions  on  (0,«),  the  right  layer  terms  as  functions  on  (-  *,0]. 


We 

assume  that  S_ ,  s 
n  p 

e  cVtO,")2  x  R2  x  [0, yXq]  )  and  that 

X  <  XQ.  Inserting  into 

(2.11), 

comparing 

0(1)-terms  and  evaluating  away  from  x  «  ±1,  X 

gives  the  reduced 

problem 

(or  order 

zero) s 

(a) 

0  =  n-p-D(x) 

(b) 

n'  “  n  ♦'♦J 

n 

(3.5) 

(c) 

p  “  -p  *'-Jp 

>  -1  <  x  <  1  . 

(d) 

■  v^p-vv0’ 

(e) 

J'  -  -S  (n,p , J  ,J  ,0) 

P  P  "  P 

We  have  to  expect  that  </,  n,  p  are  discontinuous  at  x  =  X,  therefore  (3.5)  has  to  hold 
for  the  right  hand  and  left  hand  limits  at  x  •  X.  Evaluation  close  to  X+  gives 
the  right  (zeroth  order)  (internal)  layer  problem 


(a) 

♦  » 

n-p 

(b) 

I 

n  = 

.  _  X 

(n+n(X+) )+ 

(3.6) 


-(p+p(X+>)* 


0  <  a  <  m 


i 

(d>  J 


n 


0 


(e)  J 


0 


(•••  denotes  differentaition  with  respect  to  the  corresponding  fast  variable  in  the 


sequel)  and  evaluation  close  to  X-  gives  the  left  (zeroth  order)  (internal)  layer  problem 

*  »  * 

(a)  *  ■  n-p 


(3.7) 


(b)  n  -  (n-tn(X-))* 

t  -  _  *  l 

(c)  p  -  -<p+p(x-n*  y  -  •  <  a  < 

o 


t 

(d)  j 


(e)  J  -  0 
P 


Similarly  we  obtain  the  left  (boundary)  layer  problem 

•a  _  _ 


(3.8) 


(a)  -  nt-pt 


(b)  ■  (n^+n(-1 ) ) ^ 


(c)  *  -(pjl+p(-1 ) )  ^  0  <  T  <  • 

(d)  J  -  0 


(e)  J  -  0 

P» 

The  right  (boundary)  layer  problem  is  obtained  from  (3.8)  by  substituting  n(-1),  p(  —  1 )  by 
n<1>,  p(1). 

Inserting  into  (2.19),  (3.1)  (a)  and  comparing  0 ( 1 )-coef f icients  of  1  gives  the 
matching  conditions  at  x  «  -1 


(3.9) 


(a)  n(-1)  +  nt(0) 

(b)  ?<-1)  +  pt(0) 

(c)  ^(-1)  +  *t<0) 


{ 

{ 


D(-1) 

0 

0 

-D(-1) 


*_ 


D(-1)  >  0 
D(-1)  <  0 
,  D(-1)  >  0 
,  D(-1)  <  0 


(2.18),  (3.1)(b)  gives  the  matching  conditions  at  x  ■  1 


(3.10) 


(a)  n( 1 )  ♦  n  (0 )  -  D(1) 


(b)  p< 1 >  +  pr(0)  -  0 

(c)  *(1)  +  ?r<0)  -  <>+  . 

Tha  continuity  condition*  (2.16)  give  the  interface  conditiona 


(3.11) 


(a) 

n(X-) 

+  n(0-)  - 

n(X+) 

+  n(0+ ) 

A 

(b) 

p(X-) 

+  p(0-)  - 

p(X+) 

+  p(0+) 

A 

_ 

A 

(c) 

*(X-) 

+  *(0-)  - 

*(X+) 

+  f(0+) 

t 

s 

<d) 

*(0-) 

-  *<0+) 

(e) 

J  (X- 

)  +  3  (0-) 

_  A 

»  J  (X+)  +  J  (0+) 

n 

n 

n 

n 

(f) 

J  (X- 

)  +  J  (0- ) 

__  A 

-  J  (X+)  +  J  (0+) 

p 

P 

P 

p 

From  (3.7) (d) , (e),  (3.8)(d),(e)  and  from  the  analogoua  equations  for  the  right  Jn,  J  - 
layer  terms  we  immediately  conclude  that 


(3.12) 


J  3  J  =  0,  J  S  J  =  0,  J  5  J  £0 
n  p  nt  pi  n,  P* 


since  (3.4)  has  to  hold.  No  zeroth  order  layers  occur  in  Jp.  The  current  densities 

are  the  slow  components. 

The  problem  (3.8),  (3.9)  has  been  dealt  with  in  Markowich,  Ringhofer,  et  al  (1982)  and 
it  has  been  shown  that 

(3.13,  S  "i  5  Pl  S  0 

holds.  The  same  analysis  goes  throught  for  the  right  boundary  layer  terms  and 

(3.14)  5  =  n  =  p  =  0 

r  r  r 

follows.  No  zeroth  order  boundary  layers  occur,  since  the  reduced  boundary  conditions 
for  n,  p  ((2.18),  (2.19)  with  X  ■  0)  can  be  fulfilled  by  the  reduced  solution  due  to 
<3.5)(a>. 

By  including  more  terms  in  the  expansion  (3.2)  it  turns  out  that  higher  order  boundary 
layer  terms  occur.  Similarly,  higher  order  internal  layers  occur  in  the  slow  component 

as  as  as  ^  as 

J_(J  )  if  or  t““  fr— E  or  are  not  constant  zero, 

n  p  9n  dp  '  9n  on  * 


<)i  s  n  =  p  =0 
r  r  r 
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Integrating  (3.6) (b), (c),  (3.7)(b),(c)  givaa 


(3.15) 


n(X+)(e*(0)-1)  ,  o  >  0 

A 

n(X-)(e*0>-1>  ,  g  <  0 


* 

.  f  p(X+)(a"*0,-1)  ,  o  >  0 

(3.16)  p(O)  -  < 

U<x->(a-*0,-i>  ,  a  <  0  . 

Inserting  (3.15),  (3.16)  into  ( 3. 1 1 ) (a) , (b)  gives 

(3.17)  n(X-)e*<0">  -  n(X+)e*<0+> 

(3.18)  pJX-le"^0-1  «  pU*)e~^K0*)  . 

From  (3.11)(c),  (3.51(a)  we  get  the  interface  conditions  for  the  reduced  problem 

(a)  n(X+)e^X”*”^X+*  -  n(X-) 


(3.17) 


(3.19) 


(b)  (n(X-)  -  D(X-))e^X~)”^<X+*  -  n(X+)  -  D(X+) 


and  from  ( 3. 1 1 ) (e) , (f )  and  (3.12) 


(3.19) 


(c)  J  (X+)  -  J  (X-) 
n  n 


(d)  J  (X+)  -  J  (X-)  . 

P  P 


The  boundary  conditions  follow  from  (3.9),  (3.10)  and  (3.13),  (3.14) 


f  0 

,  D(-1) 

<  0 

(a) 

n(-1)  - 

< 

(3.20) 

^  D(-1) 

,  0(-1) 

>  0 

(b) 

♦(-1)  » 

(a) 

n(  1 )  « 

D(1) 

(3.21  ) 

(b) 

*(  1 )  - 

<V  * 

Rliminating  p  from  (3.5)(b),(c)  using  (3.S)(a)  gives  the  reduce-)  equations 


(a)  *' 


D'  -  (J  +J  ) 
_ n _ p 


2n  -  D 


(b)  n*  - 


(3.22) 


(n-D)J  -  nJ  +  nD*  -1  <  x  <  X 

_  D _ P _ \ 


and 

X  <  x  <  1 


(c)  J'  -  S  (n,p, J  ,J  ,0)  X  <  x  <  1 

n  n  n  p 

(d)  J*  -  -S  (n,p, J  ,J  ,0) 

P  P  n  P 

assuming  that  D  e  C*(  [-1,X] )  rt  c'(|X,1]).  p  is  given  by  (3.51(a): 

(3.22)  (e)  p  -  n-D( x)  . 

He  obtain  the  internal  layer  problem  by  inserting  (3.15),  (3.16)  into  (3.6)(a),  (3.7) (a) 

(a)  t  »  ntX-Je*  -  p(X-)e”*  -  D(X-),  -  •  <  o  <  0 

<3-23)  S  _  ;  _ 

(b)  -  n(X+)e*  -  p(X+)e  *  -  D(X+),  0  <  a  <  • 

subject  to  the  boundary  conditions 

(c)  <i<-  •)  -  0 

(3.23) 

(d)  *<«)  -  0 

and  the  interface  conditions 

(e)  *<0+)  -  J(O-)  -  *<X->  -  <KX+) 

(3.23)  •  • 

(f)  *(0+)  -  *(0- >  . 


Because  of  (2.17)  we  require  that  the  solutions  n,  p  of  the  reduced  problem  are 
nonnegative,  that  means: 

(3.24)  Mx)  >  max( 0,D(  x) ) ,  x  e  1-1,1] 

has  to  hold.  Under  this  assumption  we  prove  a  simple  consequence  of  the  interface 
conditions  (3.19). 

Lemma  3.1.  Assume  that  (3.24)  holds  (at  least  at  X-,  X+).  Then 
D(X-)  <  0,  D(X+)  >  0  =“> 

(3.25) 

(♦(x-)  <  t(X+)  and  n(X+)  >  D(X+),  n(X-)  >  0) 


holds 


Proof.  Assume  first  that  *<X-)  -  *(X+>.  Then  <3. 19)  (a)  implies  that  n(X-t-)  -  n(X-)  and 
( 3. 19) (b)  implies  D(X-)  -  D(X+).  This  is  a  contradiction  to  the  assumption  that  D  has  a 
jump  discontinuity  at  x  -  X.  Therefore  *<X->  /  *<X+).  We  compute  n(X+)  from  (3.19) 


(3.26) 


n(X+) 


.  2<*(X-)-*(X+>) 

1  *  e 


(3.25)  follows  immediately  from  (3.26). 


□ 


We  now  give  existence  theorems  for  the  reduced  problem  and  start  with  the  simple  case 

(B). 

Theorem  3.1.  Assume  that  D(x)  >  Dg  >  0  on  f — 1,13,  D  6  C^([-1,X)),  D  e  C^([X,1J)  and 
that 

,,  S  (D(x),0,J,0,0)  3  S  (D(x),0,J,0,0)  s  0  for  all  x  e  1-1,1] 

(3.27 )  n  p 

and  all  J  e  R  . 


Then  the  reduced  problem  (3.22),  (3.19),  (3.20),  (3.21)  has  the  solution 

\ 

(a)  n(x)  2  D(x) 


(3.28) 


(b)  p(x)  3  0 


(c)  J  (x)  2  0 
P 


—  u  ,1  ds 

(d)  5  r//-i  dTTT 

T 


-1  <  x  <  X 
}  X  <  X  <  1 


u 


■1  <  X  <  X 


(e)  *(x)  = 


rx  ds 

A  .  _  .  .  0  '-1  D(s) 

—  +  *n  D(x)  -  - - : - r~  , 

U  f1  ?? 

1  J-1  D(s) 

U  f* 

-£•  +  In  D(x)  -  jj-  ~  °^)  ,  X  <  x  <  1 
U_  U_  fl  ds 

T  T  '-1  D(s) 


Proof:  Assume  ( 3.28 ) ( a) , i b) ,  (  c) .  Then  (3.27),  (3.19)(c)  imply  that  5  const  on 


[-1,1].  Fror  (3.22) (a)  we  conclude  that 
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*  ♦  t„  SisJ-  -  j  fx  -flt. 

D(-1)  n  '-1  D(s) 


-1  <  x  <  X 


♦♦ +  tn  §}#  •  7„  K  d9t  ' x  < x  « 


hold*.  Now  (3.22)  and  all  boundary  conditions  ara  fulfilled.  J  has  to  ba  calculated 

n 

fron  (3.191(a).  ((3.191(b)  la  autoMtically  fulfillsd.)  (3.281(a)  follows  than  by  using 

(3. 1 )  (c) ,  ( d )'. 

□ 

For  a  recombination  rata  R  which  dapanda  only  on  n,  p,  (3.27)  is  a  diract 
consaquenca  of  tha  equilibrium  condition.  Therefore  Theorem  3.1  holds  for  the  SRH- 
recombination  term. 

Assuming  tha  validity  of  the  asymptotic  expansions  (3.2)  (which  will  be  proven  later) 
the  theorem  implies  that  the  device  is  depleted  of  holes  (away  from  the  junction)  and  that 
the  electron  current  is  asymptotically  proportional  to  the  applied  voltage  0. 

Actually  (3.281(d)  is  a  scaled  version  of  Ohm's  law  (J^  is  the  (scaled)  resistance 

of  the  davica).  n+n  and  nn+  junctions  ara  resistors. 

Now  we  turn  to  the  case  (A).  For  simplicity  we  take  the  SRH-recombinatlon  term. 


theorem  3.2.  Assume  that  D(x)  <  0  on  l-I.X),  D(x)  >  0  on  (X,1),  |D(x)|  »  DA  on 
1-1,11  and  that  O  «  C1!)-),!]),  Dec’ux.ll).  Let  Sn,  Sp  be  given  by  (2. 14)  (SRH)  . 


Moreover  assume  that 


(3.29) 


e  <  p,  p  sufficiently  small 


holds e  Then  the  reduced  problem  (3.22),  (3*19),  (3.20),  (3.21)  has  a  locally  unique 
solution  (in  c’u-I.XlU  (X,1])4)  which  fulfills  (3.24)  and 
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(a)  *(x) 


-  tn  gj&L  +  ♦(,),  -1  <  <  x 


*+  ♦  tn  Ijfj-  ♦  tfx),  X  <  X  < 


*_-<v 

♦(x)  -  0(a  )  on  [-1,1] 


(3.30) 


+_-*+ 

0(e  )  >  0 


-1  <  x  <  X 


(b)  n(x) 


.  D( x )  +  Ole  )  ,  X  <  x  <  1 


♦-"♦a. 

(c) 

J  <x) 
n 

*  0(e 

)  ,  -1  <  x  <  1 

(d) 

J  (X) 

-  0(e 

)  ,  -1  <  x  <  1  . 

Of  course  p(x)  Is  then  given  by 


(3.30) 


(e)  p(x)  ■< 


-D(x)  +  0(e  )  ,  -1  <  x  <  X 


0(e  )  >  0 


,  X  <  x  <  1  . 

Proof;  We  introduce  (  as  a  new  dependent  variable  ( instead  of  ♦)  and  obtain  from 
( 3.22 ) (a) 


2DD 1  -  (J  +J  )D  -  2nD’ 

(a)  ♦'  -  - 2 — E - - -  ,  X  <  x  «  1 


(3.31) 


(b)  ♦*  = 


D ( 2  n-D ) 

2nD'  -  D(J  +J  ) 

_ a,  £ 

D ( 2 n-D ) 


,  -1  «  x  <  X 


Boundary  conditions  for  t  are 

(3.32)  ♦(  1 )  =  *(-1)  -  0 

and  the  interface  conditions  ( 3. 19 ) ( a ) , (b)  transform  to 


(3.33) 


.  .  *-~K  P(-1)D(1)  ♦(X-)-*(X+)  - . 

U)  0  doTodu+T  e  n(x+)  '  n(x->  =  ° 

(b)  e~  *  )u(xil  e^(X-)'^(X+)  (n(X-)  -  D(X-))  -  (MX*)  -  (XX*))  =  0  . 


* ; 


-4 


N 
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I 


We  set  w  =  e  and  z  =  (♦  ,n  ,J  ,J  ),  where 

0  _  0  0  no  p0 
^50  on  [ - 1 , 1 1 


(3.34) 


n0(x) 


0  ,  -1  <  x  <  X 

D  (  x ) ,  X  <  x  <  1 


J  (x)  =  J  (x)  =0  on  [-1,1] 


and  write  the  problem  (3.31),  (3.22)<b)  -  (d),  (3.20)(a),  (3.211(a),  (3.32),  (3.33), 

—  14 

( 3. 19 ) (c ) , ( d)  in  operator  form  F(w,z)  =  0  where  F  :  [0,»)  x  (C  ([-1,1])  ♦ 

4  8  i 

(C  ([-1,1])  x  R  ).  C  <[-l,1]>  is  the  space  of  functions  of  which  are  i-times 

=  <3>, 


continuously  differentiable  on  [  —  1 , X )  and  on  ( X , 1 ]  and  lim  f  (x)  exists  for 

x+X± 

r  ( j ) 

0  <  j  <  i.  The  space  is  equipped  with  the  norm  I  f  I  .  =  £  sup  If  (x)| 

j=0  xe[-i,i] 
xyx 

Obviously  F(0,Zg)  =  0. 

We  investigate  the  equation  D^FtO.z^y  “  (f,0)  (where  D^F(0,Zg)  denotes  the 
Fr4chet  derivative  of  F  with  respect  to  z  at  (0,zQ))  for  y  =  (y 1 ,y2 ,y3,y4 )  and 

obtain 


(3.35) 


(a)  y' 


2D* 

2 


D' 

D 


r-  o-l 


(3.35) 


(b)  y' 


2D 1 
_2 


y  +  f,  X  <  x  <  1 


0  -  =- 


y  +  f,  -1  <  x  < 


N 


i a 


1 
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j 


The  boundary  and  interface  conditions  are 


where  f  e  <CX< [-1 , 1 ]) )4. 


(c) 

-  y2(- 

1)  “ 

(d) 

y1<1)  “  <*3,  y2M) 

'  °4 

(e) 

y3(xt) 

-  y3(x-)  - 

“5 

(f) 

y4(X+) 

-  y4(x-)  - 

“6 

<g> 

y2<x-> 

‘  *7 

(h) 

y2(x+) 

“  °8 

with  a.  e  R,  i«1,...,8. 

i 

Because  of  the  Fredholm  alternative  D  F(0,z.)  is  one-to-one  and  onto  iff  it  is 

Z  0 

to-one.  Therefore  we  only  have  to  show  that  the  homogeous  problem  (3.35)  (0*0,  f  30 ) 

the  unique  solution  y  2  0. 

From  (3.35) (a)  we  get 


(3.36) 


0 


y3 


1 

0y2 


/X  <  x  <  1 


i 


Therefore 


and 

(3.37 ) (a) 


y4> 


A 


( In  D)'y^ 


By4 


0,  X  <  x  <  1 


holds.  The  boundary  conditions  are 


y4 


) 


<3.37)(b)  y^(1)  =  y^(X+)  =  0  . 

Since  6  >  0  the  maximum  principle  implies  that  y^  =  0  on  [X,1]  and  therefore 

y3  =  y2  =  y,  =  0  on  IX,  n. 


one- 

has 
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(3.35) (b)  gives 


(3.38) 


.  D' 

y2  =  '  5"  y2  +  y3 


y3  6  y2 


)  -1  <  x  <  X 


1  y2 


Again  we  get  a  second  order  problem 

(3.391(a)  y"  +  Un  |D|  ) 'y^  -  -jj  y3  =  i,  -1  <  x  <  X 

(3.391(b)  yjl-D  =  y^  (X- )  =  0 

and  the  maximum  principle  yields  y^  =  0  on  (-1,X).  y^  =  y^  =  y^  =  0  on  [  —  1 , X ] 
follows  immediately. 

Therefore  D  F(0,z_)  is  an  isomorphism  and  since  D  F(w,z„)  is  uniformly  Lipschitz 
z  0  z  u 

continuous  the  implicit  function  theorem  assures  that  there  is  a  locally  unique  solution 
z  =  z(w)  of  F(w,z)  =  0  for  w  e  [0,wQ  sufficiently  small.  Since 

F(w, zQ  )  =  0(w) 

we  get  lz(w)  -  zQ<x  1  =  O(w)  .  To  show  that  this  solution  z(w)  fulfills  (3.24)  we 
compute  the  first  order  term  z 1  of  the  expansion 

OB 

z(w)  ~  l  vi  z. 

i=0  1 


as  a  solution  of  the  equation 

DzP(0'V*1  =  •DwF(°'I0)  ‘ 


solves  (3. 35 ) (a) , (b)  with  f  =  0,  fulfills  the  interface  conditions  (setting 


i 


and  the  boundary  conditions 


n  (X-)  -  aiiUEIIi  >  o 
r  ’  d(x-) 


-  -  “s8!P  >  ” 


J  (X+)  -  J  (X-)  -  0 


J  (X*)  -  J  (X-)  -  0 


*,(-1)  -  *,(1)  -  0 

n,  (-1 )  -  n^  ( 1 )  *  0  . 

n. ,  J  ,  J  fulfill  (3.36)  and  (3.38)  and  tharafora  J  fulfills  (3.37) (a)  on  (X,1) 
1  n1  p1  P1 

subject  to  the  boundary  conditions 

jp,(U  -  o-  Vx>)  -r  'bTxTr1  < 0  • 

The  maximum  principle  implies  that  J'  ■  -  4  n.  is  negative  on  [X,1),  such  that 

P,  8  1 

Uj  >  0  on  [X, 1 )  holds. 

Similarly  we  obtain  n1  >  0  on  <-1,X).  Since  the  zeros  x  “  +1,  -  1  of  n^  are 
simple  zeros,  we  obtain  (3.24). 


The  biggest  restriction  of  the  Theorem  3.2  is  the  required  smallness  of  e  .He 


obtain  from  (3.11(a)  and  (2.24) 


(3.40) 


U  ,  UblU> 

il  U  U 

T+  t  t 

e  *  e 


Therefore  (3,29)  holds  if  there  is  a  constant  K  >  0  sufficiently  large  (but  independent 
of  A)  such  that 

u  ,UbiU)| 

(3.41)  —  <  ~ - K 

T  T 

holds.  The  applied  voltage  U  has  to  be  sufficiently  smaller  than  the  absolute  value  of 


the  built-in  voltage  ( low- inject  ion  condition). 
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Numerical  calculations  (given  in  Section  5)  demonstrate  that  (3.41)  is  not  necessary 
for  the  existence  of  a  reduced  solution. 

From  (3.30)(c)(d)  we  get  the  reduced  voltage-current  characteristic 


(3.42) 


|J  I  + 

n  p 


0  ( 


4.4 


D(1)|D(-1)| 


U 

e  T(1  +  Ot'f4!4)))  . 


In  the  case  (A)  (pn- junction)  the  total  current  density  depends  exponentially  on  the 
applied  voltage  (see  also  Sze  (1969)).  (3.42)  should  be  compared  to  the  corresponding 

result  (3.28 ) ( c) , ( d)  for  the  case  (B). 

Theorem  3.2  can  easily  be  generalized  to  more  general  functions  Sn,  Sp  which  do  not 

depend  on  J  _,  J  .  (3.30)  holds  without  change  for  S  =  S  =0. 

n  p  n  p 

Now  we  turn  to  the  internal  layer  problem  (3.23).  He  prove 
Theorem  3.3:  Case  A:  Set  D(X+)  >  0,  D(X-)  <  0.  Then,  if  the  reduced  problem  (3.22), 
(3.19),  (3.20),  (3.21)  has  a  solution  fulfilling  (3.24),  the  internal  layer  problem  (3.23) 
has  a  unique  piecewise  monotone  solution  i|/  which  fulfills: 


0  <  +<o)  <  C.  exp(  ( 1-5  )/  n(X-)  +  p(X-)  a  +  d  /  <t(x+>  '-  *(*-">) 


(a) 


<  -E./V( 


for  a  <  -Ej/  *(X+)  -  ip(X-) 


0  <  exp(  (-1+6  )/n(X+)  +  p(X+)  a  +  Dfi/  ifi(X+)  -  ip(X-)) 

(b) 

for  a  >  E J  ?(X+)  -  Y(X-) 

for  every  6  >  0  where  Cj>0,  Dj>0,Ej>0  depend  on  6  but  not  on  i(i(X±)  if 
i|i(X+)  -  1/ (X- )  is  sufficiently  large. 

Case  B:  Let  D  fulfill  the  assumptions  of  Theorem  3.1  and  let  the  reduced  solution  be 
given  by  (3.28).  Then  the  internal  layer  problem  (3.23)  has  a  unique  piecewise  monotone 
solution  which  fulfills 
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(3.44) 


(a)  lt<0)|  <  c6  exp((-1-«)/D(X-)  a),  a  <  0 


(b)  | a)  1  <  cfi  exp( (-i+«)/d<x+)  a),  a  >  o 


for  6  >  0,  where  Cg  only  depends  on  4. 

Piecewise  monotone  means  monotone  on  (-  *», 0)  and  on  (0,«). 
Proof i  For  any  piecewise  monotone  solution  f  of  (3.23) 


(3.45) 


(a)  sgn  1>(0+)  »  -sgn  *<0-),  (b)  *(0+>  t  0,  <K0-)  t  0 


has  to  hold.  This  follows  from  the  monotonicity  and  from  (3.23 ) (c) ,  ( d) ,  ( f )  since 
•  •  ^ 

♦(0+)  **  t|i(0-)  =  0  would  imply  if  =  0  (because  n(X+)  -  p(Xj)  =  D(X±)  holds)  which 

A  A 

contradicts  (3.23)(e)  because  of  (3.25)  and  (3.2S).  <i(D  +  )  =  0  (or  i(i(0-)  =  0) 
contradicts  (3.23)(e),  too. 

Only  two  possible  cases  remain: 

(I)  J(0+)  >  0,  *(0-)  <  0 

(II)  ♦( 0-)  <  0,  *(0+)  >  0  . 

A 

In  the  case  (I)  t  has  to  be  monotonically  increasing  on  (-  »,0)  and  on  (0,»),  in  the 
case  (II)  4*  is  decreasing  on  both  intervals.  In  the  case  (I)  we  derive  from  Fife  (1973, 
Lemma  2.1)  that  every  piecewise  monotone  solution  of  (3. 23 ) (a) , ( b) , (c) , ( d)  fulfills 


(  3 . 46 ) ( a ) 


^(0+)  dt 
+( a)  /2G(  t'  ) 


,  0  >  0 


(3.46) (b) 


G(t)  -  /q  ( n(X+ )es  -  p(X+)e”8  -  D(X+))ds 


n(X+)(eT-1)  +  p(X+)(e~T-1)  -  D(X+)t 


holds  ((2.23)(a)  fulfills  all  assumptions  of  Lemma  2.1  in  Fife  (1973)  because  of  (3.24)) 


( 3.47 ) (a) 


0  ,  /*<<”  _Jii_  ,  0  <  o 

4*(0- )  /2F(  X) 


( 3 .47 ) ( b) 


F ( T )  -  /J  (p(X-)es  -  n(X-)e'S  +  D(X-))ds 
»  p(X-)(eT-1)  +  n(X-)  +  D(X-)t 


(also  (2.3)(b)  fulfills  the  necessary  assumptions  after  setting  X  ■  -o). 

Differentiation  of  (3.46)(a),  (3.47)(a)  gives 

(3.48)  ii(Ot)  -  -/ 2G(Ji(0+)),  J(O-)  =  -/2F(-i(0-)) 

(3.23)(e),  (f)  gives  the  equation  for  4K0  +  ) 

(3.49)  G(Jl(0+))  =  F(<i(X-)  -  ) (X+)  -  <l(0+)) 
which  can  be  solved  (uniquely)  by  using  (3.461(b),  (3.47 ) (b)  giving 


(3. SO ) (a) 


;(0+,  ,  P(X-)(»(X+)-»(X-))  +  (n(X-)-n(X+))  +  (;(X-)-;(X+)) 
y  '  D(X-)-D(X+) 


(3.231(e)  implies 


.  50(h,  2,0  .  ,  P(X+)  (WX+l-WX-)  )  +  (n(X-)-n(X+)  )  +  (p(X-)-p(X+)) 

(  '  ¥  "  D(X+)-D(X-) 

In  the  case  (II)  we  proceed  analogously  and  obtain  the  same  formulas  for  <K0+),  iJKO-). 

A  * 

Therefore,  a  unique  piecewise  monotone  solution  of  (3.23)  exists  iff  <K0+),  <K0-)  as 
given  by  (3. SO),  have  appropriate  sign  (and  are  not  zero). 

In  the  case  ( B )  (3.50),  (3.28)  give 


(3. SI  ) 


(a)  *(0+) 


D(X- )  In  +  D(x_)  '  D(x+) 

D(X+)-D(X-) 


(b)  *<0->  = 


D(X+  )  In  |  +  D(X-)-D(X+) 

“  D (X+ ) — D ( X— )  “ 


Setting  y  =  jjj-j — *  >0  we  define 

(a)  *(0+)  =  f  (y )  =  *n*~(*~1 1 

1  y- 1 

(3.52) 

(b)  *(0-)  =  f  (y)  =  . 

2  y-  1 

Obviously  £ny  -  (y-1)  <  0  holds  for  y  >  0,  y  /  1.  Also  y£ny  -  (y-1)  >  0  for  y  >  0, 
y  /  1  holds  and  we  find  that  sqn  f^(y)  -  -r>qn  ^  ^  ^  for  Y  >  0# 
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y  f  1 .  Ha  dariva  frost  (3.52) 


(3.53) 


(a) 

♦  (0t)  >  0,  *(0-)  <  0 

iff 

D(X+) 

<  D(X-) 

(b) 

*((>♦)  <  0,  J(O-)  >  0 

iff 

D(X+) 

>  D(X-) 

In  tha  casa  (A)  va  axpraaa  n(X±),  p(X±)  in  terms  of  W**)  using  tha  lntarfaca  condition 


(3. 19)(a),  (b)  (aa  in  (3.26))  getting 
_  _  >-)-*<*♦)_, 

(n(X-)-n(X+) )  +  (p(X-)-p(X*))  -  2 -  (D(X-)  ♦  D(X+)>.  He  aat  z  -  14X-)  -  #X*> 

e*(X-)-*(Xt)+t 

and  obtain 

-  h  (z)  .  h  (*) 

(3.54)  (a)  »(0+)  ■  5(57 )"-d (77) »  <b>  »(°~)  -  5^r-p(iT) 

where 

(3.55)  (a)  (z)  -  D(X+)g2(z)  +  DtX-Jg,!*) 

(3.551(b)  h2(z)  «  D  ( X-  )  g2  ( z )  +  DtX+lg^z) 

with 

< 3 . 56 )(.-»)  g,  (z)  -  e*  -  1  -z(e*+1> 

( 3. 56 )  ( b)  g2(z)  -  e*  -  1  . 

we  restrict  to  z  <  0  since  z  »  0  cannot  occur  because  of  Lemma  3.1  g2  <  0  for  z  <  0 

and  a  simple  computation  shows  that  g^(z)  >  0  for  z  <  0.  Since  D(X+)  >  0,  D(X-)  <  0 

we  obtain  h^z)  >  0  and  h2(z)  <  0  for  z  <  0.  Therefore 

(3.57)  (a)  *(0+)  <  0,  (b)  i*0->  >  0 

follows  in  the  casa  (A). 

Now  the  existence  theorem  is  settled  in  both  cases,  the  decay  statements  (3.43), 

(3.44)  still  have  to  be  shown. 

In  the  case  (B)  the  equation  (3.231(b)  reads 

<*  -  g(i|/>  D(X+)(e*  -  1)  . 

g'(0)  -  D ( X+ )  holds  and  (3.441(b)  follows  from  Fife  (1973,  Lemma  2.1).  (3.44)(a)  is 

derived  in  the  same  way. 

In  the  case  (A)  we  have  to  keep  in  mind  that  t(>(X+)  -  i|t(X-l  can  be  large  (see  Theorem 
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3.2),  which  implies  that  <i(0+),  4»(0-)  are  large.  Therefore  we  need  estimates  which  are 
uniform  for  large  <»(X+)  -  4KX-). 

The  proof  follows  the  lines  of  the  proof  of  Lemma  2.1  in  Fife  (1973).  We  set  4“  ~4< 


in  (2.23){b),  call  the  (new)  right  hand  side  of  (41)  and  compute  /f '  (0 )  * 

/F  . 

2 


(X+)+p(x+).  It  is  easy  to  show  that 
(3.58)  f  ( s )  >  ( /f 1  ( 0  )  -  6)‘a,  0  <  s  <  6 


holds  for  6  >  0  sufficiently  small.  Since  f  is  increasing  we  get 


(3.59) 

for  Y  >  S.  Therefore 


F(T)  -  /*  f  (  s )  ds  >  S>  6  +  f  (  6 )  (  t-  6 ) 


(3.60)  (✓f  *  (0)  -  6)  r*<0+)  — —  ■  <  D./  li>(0*)| 

•  4  /F7T7  4 

holds  for  1 4(0+)  |  sufficiently  large,  where  Dj  >  0  is  independent  of  4<0+).  From  Fife 
(1973,  Lemma  2.1)  we  obtain 

(3.61)  |  4>(  a)  |  <  «  exp(-(1-«)/f'(0)c  +  Dfi/|4><0+)|) 

D«  r~- - 

for  o  > - /  1 4i(0+)  |.  (3.431(b)  follows  from  (3.50(a).  The  proof  of  (3.43)(a)  is 

/f  ’  (  0  ) 


analogous. 


□ 


Similar  estimates  holds  for  the  derivatives  of  41-  n,  p  have  to  be  computed  using 

(3.15),  (3.16).  In  the  case  (8)  p  5  0  holds  (since  p  i  0). 

J  i 

If  the  interface  condition  (3.231(f)  is  changed  to  4’(°+)  “  4,<0-)  «  0(  X)  then  the 
layer  solution  t|*  changes  at  most  by  0(  \J  4 f(X+)  -  4,(X">>  (in  the  max-norm).  This 
follows  by  applying  the  implicit  function  theorem  to  the  perturbed  equation  (3.49).  This 
will  be  needed  for  the  existence  proof  in  Section  4. 

The  width  of  the  internal  layer  at  x  =  X  can  be  computed  from  Theorem  3.3. 

In  the  case  (D)  we  obtain 
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(3.62) 


(a)  d+(X) 

(b)  d  (X) 


— 2 —  1 

1  in( — - 

— )l) 

/D(X+) 

✓d(X+) 

X 

X 

—  -  1 

tn( — - 

— )l) 

✓D(X-) 

✓d(X-) 

X  ♦  0+ 

X  ♦  0+ 


where  d+(X)(d_(X))  denotes  the  width  of  that  part  of  the  layer  which  is  right  (left)  of 
X. 

For  the  case  (A)  we  obtain 


(3.63) 


(a)  d+(X)  "Of  X~  (/*(X+)-*(X-)  +  iJtnf  —  — - 'll  )) 

/n(X+)+p(X+)  /  n(X+)+p(X+) 


(b)  d_(X)  *  of  z - (/  *(X+)-*(X-)  +  |*n(- 


:)l ))  • 


fn(X- )+p(X- ) 


n(X- )+p(X- ) 


If  the  low  injection  condition  (3.41)  holds  Theorem  3.2  gives  a  (physically  relevant) 
solution  of  the  reduced  problem  and  (3.30),  (3.63)  give 


(3.64) 


(a)  d+(X) 

(b)  d_(X) 


These  asymptotics  are  uniform  as  U  ♦ 


■».  axlstencs  of  Solution* 

In  this  section  we  prove  the  axiatenca  of  aolutiona  of  tha  singularly  perturbed 
problem  (2.11),  (2.12)  using  tha  asymptotic  expansions  (3.2). 

At  first  we  deriva  an  a  priori  aatiaate  on  tha  number  of  carrier  pairs  valid  for  the 
cases  (A)  and  (B). 

Theorem  4.1.  Every  solution  of  (2.11),  (2.12)  which  satisfies  (2.16),  (2.17)  and  for 


which  J  ,  J 
n  p 

do  not  change 

sign 

in 

[-1,11 

fulfills 

(4.1) 

U  >  0 

<—> 

J  > 

0, 

J  >  Ol 

J  )  0,  J 

I  0  on 

n 

P 

n  p 

(4.2) 

0  «  0 

<**> 

J  5 

J 

=  0 

n 

P 

(4.3) 

U  <  0 

<—> 

J  « 

0, 

J  <  Ol 

J  |  0,  J 

t  0  on 

n 

p 

n  p 

M 

1EL 

A  A 

UT 

A 

(4.4) 

***** 

Y  1  e 

T  < 

n(x)p(x)  e  Y  X 

»  T 

e  ,  x  e 

[-1,1] 

The  proof  is  completely  analogous  to  the  proof  of  Theorem  4.1  in  Markovich,  Ringhofer, 

Selberherr  and  Langer  (1982)  and  requires  only  the  equilibrium  condition  on  the  scaled 

recombination  rates  S  .  S  .  For  U  «  0  the  currant  densities  J  ,  J  vanish  and  the 
r.  p  n  p 

4  4 

device  is  in  thermal  equilibrium.  The  np-product  is  constant  Y  1  throughout  the 
device . 

The  estimate  (4.4),  the  equilibrium  condition  and  the  continuity  of  Sn,  Sp  imply 

lyl 

that  Sn,  Sp  are  small  along  a  solution  when  is  small.  In  particular,  for  tha  SRH- 

T 

recombination  rate  (given  by  (2.14)) 


(4.5) 

holds  along  every  solution  of 

(4.6) 

|0[ 

for  sufficiently  small 


JilL 

2  2  Um 

|S|  “  IS  I  «  (e  -  1> 
n  p  i  ts 

(2.11),  (2.12).  Therefore  it  is  intriguing  to  set 

S  =  S  =  0 
n  p 


We  now  give  existence  proofs  for  (2.11),  (2.12)  in  this  case. 


For  the  simple  case  (B)  we  show 


Thaorsn  4.2.  tiMN  that  D  fulfill*  the  assivptions  of  Theorem  3.1  and  that  (4.6) 


holds.  Than,  if  0  <  X  <  1Q  and  |— I  <  p  for  p  auffieiantly  saall  but  indapandant  of 

T 


X  holds,  tha  problaa  (2.11),  (2.18),  (2.19),  (2.26),  (2.27)  has  a  solution  which  fulfills 
(2.16),  (2.17)  and 


(4.7) 


*(x,X)  -  *<x)  ♦  *(~)  ♦ 


0(X) 


(4.8) 

(4.9) 

(4.10) 

(4.11) 


n(x,X)  «  d(x)  ♦  n(^j^)  ♦  O(X) 


p(x,X )  -  O(X)  (>  0) 


3  (x,X )  -  J  (x)  ♦  O(X) 
n  n 


J  (x,X)  -  0(X ) 

P 


uniformly  on  ['1,1]  whara  4,  3  ara  givan  by  (3.28)(a),  (d),  4  is  as  in  Theorem  3.3, 

n 


case  B  and  n  fulfills  (3. IS). 

Proof.  Tha  right  hand  sidas  of  (4.7)  -  (4.11)  ara  tha  sum  of  tha  raducad  solutions  as 

* 

given  in  Thsorem  3.1,  tha  layar  tarns  as  of  Theorem  3.3  (p  =  0  holds  sines  p  8  0)  and 


remainder  terms.  Ha  danota  thasa  ramalndars  by  E,,  E  ,  E  ,  E  ,  E  .  Inserting  into 

♦  n  p  Jp 


(2.11)  (with 

S  I  S  =  0),  using  (3.23)  and  (3.6) 
n  p 

(a) 

x’e;  -  e„  -  Ep  -  x2?- 

—  a  ^  a 

(b) 

E*  -  (♦♦4)*B  ♦  (n+n)Ef  ♦  E. 

n  n  f  j 

n 

(4.12) ( c) 

K  ■  -[♦+*),B  -  B  -  EE!  ♦ 

P  P  Jp  P  V 

(d) 

O 

a 

c 

tf* 

(a) 

O 

a 

i 

whara  tha  functions  4^ ,  4j  satisfy 

(4.13)  I  (^  (s,X )  |ds  -  0(X),  i  -  1,2  . 

Inserting  into  tha  boundary  conditions  (2.18),  (2.19),  (2.26),  (2.27)  shows  that  the 


boundary  values  for  E^,  En,  Ep  at  x  •  fl  ara  O(X). 


He  define  the  operators: 


(4. 13) 


where  0(x) 


(a)  (Hng)(x)  =  exp(  ^(x)  +  p(  0(x) )  -  p(s)  -  p(  o(s>  )g(s)ds 

(b)  (Hpg)(x)  -  exp(-*(x)-*(0(x))  +  i>(s)  +  i(  d(s) )  )g(s)ds 

-X 

r—  and  rewrite  (4.12)  (b),  (c)  as  integral  equations 


E  -  expHi(x)  +  *(o(x))-p<-1)-*(o(-1))E  (-D  + 
n  n 

(H  (n+n )E'  )  ( x)  +  (H  E,  )(x>  +  (H  E  E'Xx) 
n  nJ  n  n  * 

n 

+  (H  ♦.Xx) 
n  i 


(4.141(b) 


E  »  exp(-*i(x)-p(o(x))  +  *»(  a(-1  )  )E  (-1)  - 
P  P 

-  (H  Et  Xx)  -  (HE  E'Xx)  +  (H  i)(x)  . 

p  J  pp  ♦  p  2 

P 


From  (4.121(c),  (d)  we  get 
(4.151(a) 


const,  on  [-1,1 ) 


( 4 . 15 )  ( b ) 


E  3  const*  on  [-1,1] 

J 

P 


since  Ej  ,  Ej  ec([-1,1]). 
n  p 

Because  of  (4.13)  and  since  |E_(-1)|  *  0 ( X) ,  |E  (  —  1)1  =  0(1)  we  obtain  from  (4.14) 

n  p 

(a)  E  =  H  (n+n)Ei  +  E  H  1  +  H  E  E’  +  0 { X ) 
n  n  ^  J  n  n  n  p 


(4.16) 


(b)  E  -  -E  H  1  +  H  E  E'  +  0  (  X) 
P  Jp  P  P  P  ♦ 


Partial  integration  and  (3.23)  give 

(4.17)  (H  ( n+n )E '  )  ( x)  -  <n+n)E.  -  (H  J  E.Xx)  +  XG  ,EJ 

n  p  n  n  ♦  n,  X 

where  G  .  :  C([-1,1))  *  C(  [  —  1,1] )  is  uniformly  bounded  (in  X).  The  continuity  of 
n, 

n+n  at  x  =  X  was  used  for  the  derivation  of  (4.17). 

From  (4.121(b)  we  derive,  after  partial  integration 


(4.18)  H  E  E'  =  E  E  -  H  nE'E,  -  H  E  E  -  H  E  E'E,  +  XF  ,E\E  ,  +  0  (  X) 

n  n  f  n  i|>  n  ip  ip  n  J  ip  rt  n  n,  X  ip  ip 

n 
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where  F 


t  C( 1-1,1])  ♦  C( 1-1,11)  is  uniformly  bounded.  Another  pertiel  integration 


(4.19) 


HnnEjE^  -  j  nE^  -  •*/“,  (e  *n) 'E*(s)ds  + 


Combining  (4.16)  -  (4.19)  gives 


(4.20) (a) 


En  -  («♦«>**  -  ♦  Bj  Hn1  ♦  x(E#,E;.En.E  > 


and  proceeding  analogously  for  Ep 


(4.20) (b) 


E  “  -E  H  1  +  A  \ (E . ,E!,E  ,E_  ) 
P  PP  P»**VPJp 


where  A  .  ,  A  .  are  nonlinear  operators  from  (C([-1,1]))4  into  C((-1,H)  which 
n.X  p.A 


fulfill 


(4.21) (a) 


(4.21 ) (b) 


'V'VVVV'm.11  4  C1<0(X)  +  'V (-1.13 (X+X,V  [-1.1] 
+ 1 V  i-i#i)  +  'Vi-1.11  +  'V  t-i.u'V  1-1,11 1 } 

,DXq.X  ‘  C2<X  +  *V  [-Ml 

♦'Vl-I.U  +,BJq'M.1] 

+  'Vm.IJ'Vm.U  +  *Bq'  I— 1.1]  '*♦'  [-1,1] 

♦  '*+' i-i.u'V  (-Mi  +  x,v  1-1.11* 


where  C,,  Cj  are  independent  of  X,  "D"  denotes  the  Frfechet  derivative  and  q  -  n.p. 
The  constants  Jn,  Jp  can  be  determined  from  the  boundary  conditions  for  En,  Ep  at 


x  -  +1« 


(4.22) (a) 


(4.22) (b) 


‘VnV<,)  *  An.X%'E*'En'EJn,,1) 

(H  1 )  ( 1 ) 
n 

Ap,X  ^E4i'EiJ'En'EJn^ ( ' » 

(H  1 ) ( 1 ) 

P 


We  remark  that  (Hn1)(1),  (Hp1)(1)  are  bounded  away  from  rero  uniformly  in  X. 

We  regard  E^,  En,  Ep,  Ej  ,  Ej  as  dwelling  the  space  Ajj  -  c’u-I.U)  D 
C*([-1,1]>  which  is  equipped  with  the  norm  I  f  I  x  -  If!  ,,  ,,  +  x2,f"1  j.,  x)g(X,1]' 


implies  that  the  mappings  defined  by  the  right  hand  sides  of  (4.20)  are  contractions  in 
spheres  of  radius  0(A)  (centered  at  0)  and  therefore  (4.20)(a)t  (b)  can  be  resolved 
with  respect  to  En,  Ep  reap.: 


(4.231(a) 


(4.231(b) 


(H  1 ) (x) 

<n+n)%-  H„W  (VnV("  oTiTdT 

n 

EP  =  Qp.X(E*'E*)  • 


The  operators  0  , ,  0  .  fulfill  the  estimates  (4.21)  when  IE  I.  .  ,,,  IE  1.  . 

n, a  p,A  n  (-1,1)  p  (-1,1) 

lEj  *  ^ | >  lEj  I  ^  j  are  substituted  by  the  radius  of  the  sphere  in  which  the 

n  p 

contraction  mapping  theorem  is  applied. 

Inserting  (4.23)  into  (4.121(a)  gives 

,  _  (H  1 )  ( x) 

A  E”  -  (n+n )E  -  (H  J  E.)(1)  7—":  .  -  (M  J  E.)(x) 

♦  ♦  n  n  ♦  (H  1  )  ( 1 )  n  n  ♦ 


(4.24) 


+  “n^W  '  nP,A(VEV 


■  X  (T,  -1  <  x  <  X  and  X  <  x  <  1 


subject  to  0 ( X )  boundary  conditions  for  E  at  x  =  tl . 

V 


Since  n+n  is  positive  and  continuous  on  [-1,1]  the  boundary  value  problem 


(4.25) 


X  y"  -  (n+n )y  =  f(x),  y ( —  1  )  =  y_,  y(1)  »  y+ 


has  a  unique  solution  y  e  Ax  for  all  y_,  y+  e  R,  f  e  C  ([-1,1])  which  fulfills 

(4.26)  lyl^  <  const ( If I |_g +  +  ly+l*  * 

Since  all  estimates  so  far  are  uniform  for  I"- I  <  p  and  since 

T 

(4.27)  IH  J  E.l.  ,  ,,  <  const  plE.I.  ,  ,, 

n  n  [-1,1]  ♦1-1,1] 

holds  with  const  independent  of  X  and  p  e  [  —  pQ , pQ ]  (see  (3.2H)(d))  the  contraction 
mapping  theorem  with  p  sufficiently  small  assures  the  existence  of  a  locally  unique  E 


V[-1,U  '0U)  • 


T  1 


The  existence  statement  of  Theorem  4.2  follows  immediately.  The  positivity  of  p  ”  Ep 
follows  by  investigating  the  higher  order  terms  of  the  expansions  similar  to  the  proof  of 
Theorem  3.2. 


This  proof  does  not  carry  over  to  the  case  A  since  then  Hn,  Hp  are  not  uniformly 

bounded  (in  X)  anymore.  At  first  we  rewrite  (2.11)  (with  S  2  s  SO)  as  a  second 

n  p 

order  problem. 


From  (2.11)(d)(  (e)  we  get 


J  2  const,  J  =  const  on  [—1,1] 
n  p 


(2.11 )(b),(c)  give 


(4.30) 


(a)  n  -  n(1)e*(x)-*(1)  +  j  e*(x>  /*  e'^’ds 

n  '  1 

(b)  p  -  p(1)e*(1,-*(x)  -  J  e-*(x)  f;  e*B,ds 


Jn,  Jp  have  to  be  determined  from  the  boundary  conditions  for  n,  p  at  x  *  -1 


(4.31)(a) 


(4.30(b) 


He  uj, 

2  2/  °T  UT, 

aii.,  *  ,_  J 

e-*,8,ds 

\  !k 

2  2  /  °T  V 
=  r  *  i*  -  ■»  ] 

f  1  1>(s> 

J_  e  ds 


(4.30) ,  (4.31)  immediately  give  n  >  0,  p  >  0, 

Without  loss  of  generality  we  set  “  -Uc  “  j  and  obtain  by  inserting  (4.30), 

(4.31)  into  (2.10(a) 


(4.32) 


X2*"  =  2y2  X2  (sinh( 


fx 

T  T  J_1e  ds 


f*  t's). 

-*  h*  ds  , 

*  p  — — — : — ))  -  o(x) 

/■I  *(s).  ’ 


for  x  e  (-1,X)  (J  ( X , 1 ]  subject  to  the  boundary  conditions  (2.20),  (2.21)  and  the 


interface  conditions  (2.16). 

At  first  we  derive  an  asymptotic  representation  of  the  voltage  current  characteristic. 

Lemma  4.1.  Let  S  =  S  =0  hold  and  assume  that  U  fulfills  the  low  imection  condition 
-  n  p 

(3.41),  that  D  fulfills  the  assumptions  of  Theorem  3.2  (Case  (A))  and 
(4.331(a)  i>(x,  X)  =  *<x)  +  *(^)  +  0<X) 

where  the  reduced  solution  +  fulfills  (3.30)(a)  and  i(i  is  given  by  Theorem  3.3.  Then 
(3.31)  implies 


(4.33) (b) 


J  = 


u_ 

4.4,  UT 
Y  X  (e  -  1  ) 


/^1 |D(s) Ids 


( 1  +  0<  |  8(  X)|  +  dx(  X)  + 


4  4 

d  (X)  +  Y  X  e  )  ) 


(4.331(c) 


U_ 

4  4 

=  Y  X  (e  -  1) 
/'  D(s)ds 


(  1  +  0(  |  6(  X)1  +  d  (X)  +  d  (X)  + 


U 

4  ,4  T, 
Y  X  e 


Proof .  We  obtain 


fl  -i|/( s ,  X )  _  fx  D(s) 

°  ' -\e  ds  "  J-i  dcT> 


.s-x 

rx  o\s)  ,  (X  D(s)  V  X  0(|  0(  X)  |  ) ds 

ds  -  ET^TT  ds  +  •'-I  e 


D(-1  ) 


+  fx  DC)  .  ...  +  A  D ( 1  ) 

+  J  1  — 77  <e  “  1  )ds  +  e  j  ~ — r  ds 

‘'“I  D(-1  )  Jx  D(  s  ) 


•  j  i  /  s-x . 

,  fi  Dd)  -*,-r> 

+  e  L  dTTF  (e  -  1)ds 


—  W/  —  - ) 

+  e  "  *  /'  £7—  e  x  0  (  (  9(  X)  |  )ds 


'x  D ( s ) 


Using  the  estimates  (3.43)  and  the  layer  widths  given  by  (3.64)  we  derive 


f\  e‘*(s'i,d3  =  £1  sfriy  ds  +  °< i °<^>»  +  o(dj»> 


+  0(y4X4e  T  )  +  0(d+(A))  . 


Therefore 


*.  2  2  2UT 
b  r  A  (e  -  e 


e  -  /l,  e-*(s'X1ds 


■*-  * 


Y4A4(eUT  -  1) _ 

e-*(s'X)dS(i  +  0(t4x4)) 


4.4.  ^  ,  V 

1  A  (e  -  1 ) 

J*j  |D(s) I ds 


(1  +0(|  0(  X) ! 


+  d  <  X)  +  d  (  X)  +  y4X4e  T)  ) 


follows.  The  proof  for  the  asymptotic  representation  of  is  analogous. 

If  D  *  D  holds  (which  implies  that  J J_  have  the  same  scaling  factors)  and  if 
n  p  n  p 

a  solution  of  (4.32)  subject  to  the  boundary  conditions  <i(-1,A)  «  <i_,  (K 1 ,  X)  *  exists 
for  which  8(A)  +  0  as  X  ♦  0+,  then  the  total  voltage  current  characteristic  of  the  pn- 
function  is  given  by 


(4.331(d)  J  £  J  +  J  S  y4A4( - - -  +  — - - - )(e  T  -  1)M  +  o(  1 ) ) 

n  P  | D( s ) Ids  J^D ( s ) ds 


The  same  asymptotic  form  of  J  can  be  found  in  Sze  (1969)  (and  other  standard  books  on 
semiconductor  physics),  however  the  derivation  used  there  heavily  relies  on  physical 


arguments. 


We  now  prove  an  existence  Theorem  for  the  case  (A)  under  a  slightly  sharper  assumption 


on  l—l  then  that  used  in  Theorem  4. a* 

T  2 
Theorem  4,3.  Let  D  fulfill  the  assumption  of  Theorem  3.2  and  D  e  C^l (-1,11  I.  Assume 


U  ..  2.2. 

T  |  tn  Y  X  I 


holds  for  some  p  >  0  sufficiently  small  hut  independent  of  X.  Then  there  is  a  solution 


ifdx.A)  of  (4.32)  subject  to  (2.20),  (2.21),  (2.16)  and  i(i  fulfills 


9 


(4.34) 


4<x,X)  “  *(x)  +  i»(^)  +  0(d+(  X)  +  d_ ( X)  +  (Y4XV) 
for  some  y  >  0  where  4  is  the  reduced  solution  as  given  by  Theorem  3.2  (with 
Sn  5  Sp  -  ^ 1  and  ♦  is  the  layer  solution  given  by  Theorem  3.3. 

Proof.  We  define 


=  In 


4  =  in 


D(1)  +  i/D(  1  )2  +  4Y4!4 


2Y2X2 


L-D(-i) 


♦  /d(^'2  ■  *-4-4 


-1)  +  4y  X  j 


Since  2y2X2  “  — — 1 ^  1 -  holds,  the  problem  (4.32),  (2.20),  (2.21)  with  U  »  0 

sinh  4+  sinh  4 

can  be  rewritten  as 


(4.351(a) 


(4.35)(b) 


(4.351(c) 


(4.351(d) 


sinh  <|» 

X  4"  =  0(1)  - j  -  D(x>,  X  <  x  <  1 

sinh  4+ 

2  Sinh  ♦o 

X  4J  =  D(-1)  - -  -  D(x),  -1  <  x  <  X 

sinh  4 

♦Q(1,X)  ■  4°,  4Q(-1,X)  =  4° 


e  c\  [-1,1]  )  . 


We  now  regard 


parameters  independent  of  X  (as  in  Section  3). 


Then  the  reduced  solution  4g  has  the  asymptotic  form  given  by  Theorem  3.2. 


(4.36)  4»q  <  x) 


■  D(x) 


D(x) 


area  sinh(  —  sinh  4+ )  =  4+  +  in  ■■■--  +  0(e  ),  X  <  x  <  1 


.  ,  ,D(x)  .  .  .0.  .0  .  D(x) 

area  srnh(  .  ;  sinh  ip  )  =  4*  “in  ~ — ~  +  0(e 


-24 


D  (  - 1  ) 


D  (  - 1  ) 


),  -1  <  x  <  X 


and  the  internal  layer  solution 


in  Theorem  3.3  but  subject  to  the*  changed 
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1  !  _ 

interface  condition  ^<0+)  -  0— )  •  X( f (X- )  f'  (Xt) ) . 

In  order  to  Investigate  stability  we  substitute  uQ 


*0 


in  (4.32).  (2.20),  (2.21) 


with  U  •  0  and  obtain  the  problem 

2  D(  1 )  **n*'  K%  1 

(4.37 )  (a)  X  u“  -  - jp  -  —  D(x)  -1  «x<X,  X<x<1 

sinh  *+ 

( 4. 37 )  ( b)  “-(-I )  =  —  ,  u  (1)  -  1 

K 


(4.37)  (c)  u  6  c’d-I.D  )  • 

He  denote  (4.37)  by  F0(u0,X,<i°,<i+),  where  F..(  *,  X,  )  t  Ax  ♦  Cx(  [-1 , 1  ]  )  x  R2  and 

investigate  the  equation 

(4.38)  LX,0V  S  DuF0(U0+U0,X'^'*-)v  ”  tf,O,0) 


where  uQ 

(4.38) (a) 


*0  *  *0 

Uq  *  — .  (4.38)  is  equivalent  to 

if  If 

%  T+ 


X  v"  -  D(1) 


sinh(*0+*0) 

sinh 


v  »  f(x). 


<x<X,  X<xC1 


(4.391(b)  v(-1)  -  o,  v(1)  -  B 

(4.391(c)  v  e  C*(  1-1,11  )  . 

-  *  .  1 

We  remark  that  +  4iq(^)  *  C 

The  maximum  principle  immediately  implies  uniqueness  of  the  solution  of  (4.39)  and  the 
Kredholns  alternative  gives  existence.  To  get  a  bound  for  the  inverse  of  (4.39) 
we  construct  the  barrier  function: 


(4.40) 

where  the  constants 


v  (x.X)  =  K  +  exp(-  P;-‘*L5") 


K  >  0,  B  >  0  will  be  determined  thereafter. 
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We  compute 

(4.4D  <L  v  )(x)  =  (-§*-5  (2f-4  -  ’Wf"  y'o'V)  - 

'  x' ‘if  -i  ) 


cosh  (</  +♦.  ) 

D(1)  - 5 —  (k  +  exp(- 

sinh  <i+ 


P(x-X) _ \\ 

2  0  O''  ' 

X 


^  (  ( Bd-x)2 


K  ♦  exp(-  K  ♦  exP(-  ^fV))  • 

x2<i£-0  x2(<.“-ij._) 


denote  the  first  component  of  (L^  ^v^Mx)  ^X  gv^)(x)  an<^  9et 

1  %/0“0  . /o  .0 

(,)  g  -  J  X/V*-  X/V*- 

(L  „v.  )(x)  <  -  -r2-^  e  ,  x  e  X - ,  X  +  - — 

X'°  b  *°-*°  2/i  2/6 


,  26  ,28<x-X)  /  6<x-x>  .  46e  _  ,  . 

(~TT  (  2  o  0  ~  deXPl~  2  0  0  "  ‘  00  '  X  e  ["1'11 

where  (f)+  denotes  the  positive  part  of  the  function  f. 

Now  we  use  the  estimates  for  given  in  Theorem  3.3  for  fixed  0  <  6  <  1.  At  first 
we  estimate 


cosh(*0+<>0) 
sinh  ^ 


D(x)e  ,  x  C-[X,1] 


1 D( x ) | e  ,  x  C- [-1 ,X) 


Since  is  monotonically  increasing  on  (-1,X)  and  on  (X,1)  this  yields 


DU)  exP«.0(-^)).  *  e  |*  *  2/j 


|D<x>|exp(-*,( - X"C 
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-1.x  -  - 


v-x  a 

(note  that  a  *  1“)*  We  no"  choose  S 


such  that 


holds  and 


Then 


/6  =  5  ■*»(»“• 


1  (1-6)/ n(X+)+n(X+)  (1-6)  /  n(X-)+p(X~) 


2D  . 
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d(x)  ,  x  e 


vA°-<? 

X  +  — “ ,  1 
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ul’n^JU)  <  -  Ke’X 
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|d(x)|  ,  x  e 
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.0 

V*- 


e  1  min  |D(x)] 

xe[-1.1] 


(L^'v.  )(x)  <  -  min (l ,  ^  r)  .  x  «  1-1. H 


holds.  The  maximum  principle  implies  the  estimate 


X,°  CX(C-1,1])XR2*C([-1,1]> 


<  const .  t  ) 


(4.39)(a)  gives 


(4.42) 


IL.1  I  <  const.  (<t-*°) 

cx([~1'1])KR 


The  Frfichet-derivative  DuFo  is  ^oca^y  Lipschitz  continuous: 

ID  F  -  D  FQ(w  ,  X,  *°,*°)l 

J  U  U  "  Ax«cx(  1-1,1)  )«2 

(4.43)  0 

const.  4,  Iw, -w,  I. 

+  1  2  A„ 


for  w, ,  w,  in  a  sphere  centered  at  u  +u.  = 


T0  T)  const 

— ; —  with  radius  — r — . 


Now  we  rewrite  (4.32)  as 


(4.44) 


.2  „  D( 1 ) 
x  u«  «  __ 


sinh(  u+  rTT“)  .0  . 

m  *  2ut  ,  Dm  .  fc<JL .  H(»+u) 

0  00  Slnh  20  0 
sinh  <i+  <>+  T  sin  h 


,  X  e  [-1,X)0(X,1] 


where  H(<i)  denotes  the  integral  operator  on  the  right  hand  side  of  (4.32).  We  denote 

*_  K  o 

(4.44)  subject  to  the  boundary  conditions  u(-1)  =  — ,  u(1)  “  —  (with  ♦  + 

*+  *+ 

♦+  *  ■  j)  by  Fy(u,X,^+,t_)  =  0.  The  FrSchet  derivative  Lx  o  ” 


°uFU  <  '•'o  +  ’('o '  X ^ *  is  <3iven  by 


D ( 1  )  sinh (  4>0+<|q+  — ) 

(L  v )  ( x )  =  (X2v" - - - — 

sinh 


D(1)  - (D_h(*„+*„)v>(x),  v(-1  ) ,  v(  1  ) )  . 


D|  (N 


where 


,DuH(V‘Vcx([-1.1))  *CX<{-1,1]>  *  con#t*  »inh  ♦♦  • 


Therefore,  if 

(4.46)  Ijj-I  <  Qfi—>  P  euff.  null 

T  ♦+"*_ 

holds,  L.V  fulfills  the  estimate  (4.42).  Moreover  D  F„  is  locally  Lipschitz 
A, U  »  “ 

continuous  and  fulfills  (4.43). 

Since  -  (♦„“♦„ ) I .  ,  ,,  <  const  |U|  a  simple  perturbation  argument  shows  that 

also  )  1  fulfills  the  stability  estimate  (4.42)  and  that  f>uFu  is  locally 


Lipschitz  continuous  around  ((4.43)  holds)  if  U  is  restricted  by  (4.46). 

.  hi 

♦: 


We  now  insert  u+u  *  into  in  order  to  compute  the  'local  error'#  We  obtain 


.•*°I  -  (4  •  .  4  -  |Wll 


fu<u+u,A,*+,*_)  -  <—  <i  +  —  ♦< o)  -  (• 


sinh(  5^-) 
_ T 


sinh 


sill  Sinh(_s_)  juhh .  aisi,,  til,  ^1±1) 


with  a  =*  Obviously 


^  \  <?.*** r,  *  v 

T  sinh  <i+ 

~  -  * 
where  J  ,  J  are  given  by  (4.31)  (a),  (b)  when  ♦  is  substituted  by  <>*■♦. 
n  p 


We  rewrite 


A  A  _  A  _  A  Sinht.y  )  _  A 

J  #'♦■*  f*  e<**ds  =  J  e'^*  /*,  e ** - £-  e’^* 

P  1  P  -1  sinh 


Since  is  positive  on  (“1#X)  and  negative  on  [X,tJ  and  monotonically  decreasing  on 

[-1,X)  and  on  (X,U  we  get 


J1  e  vda|  <  const 


a,  < 


x  e  1-1,1] 


const 


Lemma  (4.1)  with  9(A)  -  0  gives 
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U_ 
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if  the  low  injection  condition  holds. 


We  calculate,  using 


D(1) 


-  =  D(-1  ) 
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Using  Theorem  3.2  gives 
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- 2 - -  -  n(X+)  +  (D(x)  -  D(X+) )  +  0(y4X4(1  +  e  T)) 

2  sinh 


r  sl 

.  . . . . r-  -  p(X+)  +  0 ( Y4 X4 ( 1  +  e  T)) 

2D(x)sinh 


x  P  (X,1) 


2  *  r  j- 
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2D( x) sinh  iji 


x  e  t-i.x) 


2  sinh 


-  p(X-)  +  (D(X-)  -  D(  x) )  +  0(Y4X4(1  +  e  T)) 


Since  f  solves  the  layer-equations  (3.23 ) (a), (b)  we  get  for  (4.47) 


(4.48) 


U_  U_ 

•  u  “  u  * 

(Dfxj-DIX+Jle*-!)  +  0(y4X4(1  +  e  T))e*  +  0 { Y4 X4 ( 1  +  e  T))e-*  x  e  <X,1J 


(D(x)-D(X-))(e-*-1)  +  0 ( y4 X4 <  1  +  e  T))e* 
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u  * 

+  0(Y4X4(1  +  e  T )  )e”^,  x  e  f-t,X) 


♦  is  negative  on  (X,1]  and  positive  on  t-1,X)  and  since  D  e  C  ([-1,1])  we  obtain 

A 

A 

|<D(x)-D(X+))(e*-1)|  <  const.  d+(X),  x  e  (X ,  1  ] 

A 

|(D(x)-D(X-))!e‘*-1)|  <  const,  d  (X),  xc  [- 1  ,X )  . 


Also,  (3.501(a),  (b)  imply 


-45- 


e  "  «  .  —  <  const.  exp(D(»:;:D;i.- 


/  |D(X-)|  (J  ,.4,4.D(X+)-D(X-)  _  . 

*  const-  exPW)-bu-'>  iTJ<Y  x  1  •  x  e  (x-15 
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and  analogously 
*  jt—X 


T  X  .  *(0-)  .  ,  -D(X+)  U  4 .4. D(X+ )-D(X- )  _  . 

e  <  e  <  const.  «p[D(x+);0  J.  ~ J(Y  X  )  ,  x-C  I-1,X)  . 

T 

U  4  4  Y 

Therefore  the  expansion  (4.48)  is  bounded  by  const  ( d_  ( X)  +  d+(X)  +  exp(et|~|)(y  X  )  ) 

UT 

where  a  >  0,  y  >  0  holds. 

—  2  2 

These  estimates  and  e  C  <  [—1,1]  1  for  D  e  Cx  (  [  —  1,1]  )  imply  that 


(4.49) 


lF0(u+ufX,*+,*.)lc^([_nin  <  const 


T  4  4  Y 

d_(X)  +  d+(X)  +  e  (Y  X*)T 


The  stability  estimate  (4.46)  (which  holds  for  y  if  U  fulfills  (4.46)), 
Lipschitz  continuity  of  DuFy  and  (4.49)  make  it  possible  to  apply  the  version  of  the 
implicit  function  theorem  given  by  Spijker  (1972),  which  implies  that  Fyfu,  X,  t(? )  =  0 

*  X 

has  a  solution  u  which  is  unique  in  a  sphere  in  Ax  with  radius  — — —  for 

<v 

sufficiently  small  x  centered  at  u+u  and  the  estimate 

lu-u  I  <  const,  (d  <X)  +  d  ( X)  +  (Y*X4)  ) 

AX  + 

holds  for  |jy-|  <  where  p  is  sufficiently  small  but  independent  of  X. 

T  *+ 

□ 

We  remark  that  the  size  reduction  on  |jj— (  comes  from  the  interpretation  of  (4.32)  as 

T 

perturbation  of  the  equilibrium  problem  (U=0)  which  was  heavily  used  for  the  stability 
proof.  The  numerical  results  demonstrated  in  the  next  Section  indicate  that  existence  and 
validity  of  the  asymptotic  expansions  hold  under  much  weaker  restrictions  on  |;~| . 


5.  numerical  Experiments 

We  demonstrate  numerical  results  for  two  pn- junctions  in  the  high  injection  case,  that 

is  U  >  li.(X).  The  existence  Theorem  3.2  for  the  reduced  problem  does  not  hold  if  the  low 
~  bi 

injection  condition  (3.41)  is  neglected.  Both  functions  we  investigate  have  characteristic 
length  2t  ■  5  x  10  3cm,  the  doping  profile  of  the  first  pn-junction  (called  junction  I 
in  the  sequel)  is 


-0.5  x  1017cm3  , 


,„17  3 

10  cm 


,  f  <  *  <  t 


and  for  the  second  junction  (called  junction  11  in  the  sequel) 


.„15  3 

■10  cm 


,n17  3 

10  cm 


-t  <  z  <  f 


,  |  <  i  <  t  . 


Accurately  speaking,  both  devices  are  pn+  junctions.  In  both  cases  C  “  10 1  cm3.  This 
and  the  numerical  values  for  the  parameters  from  Table  1  givea  for  both  devices  using  the 


formulae  (2.6),  (2.9)j 


For  junction  I  we  obtain 


X2  -  0.4  x  10‘6  ,  y2  -  0.25 


min|C(z)  | 

z _ _  _1_ 

max|C(z)|  2 


and  for  junction  II 


min|C(z) | 

z _ 

max | C( z ) | 


For  both  cases  the  singular  perturbation  approach  seems  applicable  because  (2.22),  (2.23) 


'holds'  (the  order  of  magnitude  of  X  and  min |C( z) I /max |C( z) I  as  given  by  (5.3)  and 


(5.4),  (5.5)  resp.  are  clearly  different  and  n^ 


1010  while  C  »  101?). 


-47- 


The  built-in-voltage  (calculated  ualng  (2.25))  for  junction  I  la 

(5.6)  Ubi  -  -0.79V 
and  for  junction  II 

(5.7)  Ubi  -  -0.69V  . 

All  calculation*  described  in  the  sequel  were  performed  on  the  CDC-Cyber  74  computer  of  the 
Technical  University  of  Vienna  with  the  boundary-and- inter face-problem  solver  PASVA4 
written  by  M.  Lent ini  and  V.  Pereyra.  The  SHH-recombination  term  was  used.  Figures  1-3 
show  the  reduced  solutions  of  a  typical  high  injection  case  for  junction  I  with 
U  ■  1.39V.  The  majority  carrier  densities  (n  on  the  n-slde,  which  is  the  interval  (“,  1] 
and  p  on  the  p-side  which  is  the  interval  t-1,j))  are  larger  than  the  doping  |D| 
(except  at  the  boundaries  x  -  ±1). 

The  reduced  solutions  for  a  high-injection  case  (U  -  0.99V)  for  junction  II  are 
shown  in  Figures  4-6. 

Since  in  both  cases  the  applied  voltages  are  significantly  larger  than  the  absolute 
value  of  the  built-in-voltage,  the  existence  Theorem  3.2  for  the  reduced  problem  cannot  be 
applied.  However  the  presented  numerical  results  give  a  strong  indication  for  the 
existence  of  reduced  solutions  even  in  the  high  injection  case. 

Figure  7-9  and  10-12  show  the  solutions  of  the  singularly  perturbed  (full)  problem 
with  U  “  |Ufeil  for  junction  I  and  II  respectively.  The  internal  layer  in  the  components 
4*,  n,  p  is  clearly  visible.  The  solutions  of  the  corresponding  reduced  problems  (whose 
existence  is  also  not  covered  by  Theorem  3.2)  were  also  computed  and  they  agreed  up  to 
graphical  accuracy  with  the  full  solutions  away  from  the  layer  (see  Figure  13,  which  shows 
the  reduced  solutions  n,  p  for  function  I).  In  fact,  the  reduced  solutions  were  used  as 
starting  guesses  for  the  numerical  method  to  compute  the  full  solutions  and  convergence  was 
achieved  in  a  few  steps. 

This  indicates  that  the  asymptotic  expansions  are  valid  for  a  much  larger  range  of 
0  values  than  given  in  Theorem  4.2. 
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